1 HOK comparison of CASAL and Casal2 model configurations

This document compares the results of at least 2 CASAL model configurations (base and at least one sensitivity) and up to 8 Casal2 model configurations (3 BetaDiff, 2 CppAD, and 3 ADOL-C).

The CASAL model sensitivity 1 has a smaller tolerance value than the CASAL base model (1e-9 vs. 1e-6).

The Casal2 ADOL-C and BetaDiff low tolerance models have a smaller tolerance value than the CASAL base model (1e-9 vs. 1e-6). The Casal2 CppAD models have a tolerance value of 1e-9.

1.1 R environment

## [1] "Tue Feb 11 15:25:10 2020"
## R version 3.6.0 (2019-04-26)
## Platform: x86_64-redhat-linux-gnu (64-bit)
## Running under: CentOS Linux 7 (Core)
## 
## Matrix products: default
## BLAS/LAPACK: /usr/lib64/R/lib/libRblas.so
## 
## locale:
##  [1] LC_CTYPE=en_NZ.UTF-8       LC_NUMERIC=C              
##  [3] LC_TIME=en_NZ.UTF-8        LC_COLLATE=en_NZ.UTF-8    
##  [5] LC_MONETARY=en_NZ.UTF-8    LC_MESSAGES=en_NZ.UTF-8   
##  [7] LC_PAPER=en_GB.UTF_8       LC_NAME=C                 
##  [9] LC_ADDRESS=C               LC_TELEPHONE=C            
## [11] LC_MEASUREMENT=en_NZ.UTF-8 LC_IDENTIFICATION=C       
## 
## attached base packages:
## [1] stats     graphics  grDevices utils     datasets  methods   base     
## 
## other attached packages:
## [1] casal2_1.0     casal_2.30     devtools_2.2.1 usethis_1.5.1  rlist_0.4.6.1 
## [6] ggplot2_3.2.1 
## 
## loaded via a namespace (and not attached):
##  [1] Rcpp_1.0.3        pillar_1.4.3      compiler_3.6.0    remotes_2.1.0    
##  [5] prettyunits_1.1.1 tools_3.6.0       testthat_2.3.1    pkgload_1.0.2    
##  [9] pkgbuild_1.0.6    digest_0.6.23     memoise_1.1.0     evaluate_0.14    
## [13] lifecycle_0.1.0   tibble_2.1.3      gtable_0.3.0      pkgconfig_2.0.3  
## [17] rlang_0.4.4       cli_2.0.1         yaml_2.2.1        xfun_0.12        
## [21] withr_2.1.2       dplyr_0.8.4       stringr_1.4.0     knitr_1.28       
## [25] desc_1.2.0        fs_1.3.1          rprojroot_1.3-2   grid_3.6.0       
## [29] tidyselect_0.2.5  glue_1.3.1        data.table_1.12.8 R6_2.4.1         
## [33] processx_3.4.2    fansi_0.4.1       rmarkdown_2.1     sessioninfo_1.1.1
## [37] callr_3.4.1       purrr_0.3.3       magrittr_1.5      backports_1.1.5  
## [41] ps_1.3.0          ellipsis_0.3.0    scales_1.1.0      htmltools_0.4.0  
## [45] assertthat_0.2.1  colorspace_1.4-1  stringi_1.4.5     lazyeval_0.2.2   
## [49] munsell_0.5.0     crayon_1.3.4

1.2 CASAL and Casal2 model output

# read in CASAL MPD results
cas_mpd_filename <- file.path(base_dir, 'CASAL', mpd_run)
cas_mpd  <- casal::extract.mpd(cas_mpd_filename)
cas_corr <- casal::extract.correlation.matrix(cas_mpd_filename)
## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion
cas_conv <- get_convergence_information(cas_mpd_filename)
c1_quant <- cas_mpd$quantities

cas_mpd_sens1_filename <- file.path(base_dir, 'CASAL_sens1', mpd_run)
cas_mpd_sens1  <- casal::extract.mpd(cas_mpd_sens1_filename)
cas_sens1_corr <- casal::extract.correlation.matrix(cas_mpd_sens1_filename)
## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion

## Warning in casal.string.to.vector.of.numbers(lines[i]): NAs introduced by
## coercion
cas_sens1_conv <- get_convergence_information(cas_mpd_sens1_filename)
c1_sens1_quant <- cas_mpd_sens1$quantities
C2_subdir <- c('betadiff_casal_flags_on',
               'betadiff_casal_flags_off',
               'betadiff_casal_flags_on_low_tol',
               'cppad_casal_flags_on',
               'cppad_casal_flags_off',
               'adolc_casal_flags_on',
               'adolc_casal_flags_off',
               'adolc_casal_flags_on_low_tol')


# find the non-zero-sized parameter estimate files for the C2 models
C2_subdir <- unlist(sapply(C2_subdir, function(c_dir) if ( file.access(file.path(base_dir, 'Casal2', c_dir, params_est), 4) == 0 &
                                                           file.size(file.path(base_dir, 'Casal2', c_dir, params_est)) > 0 )
                                                         { return (c_dir) }, USE.NAMES=FALSE ))


num_C2_models <- length(C2_subdir)

C2_color <- c('blue', 'green3', 'red', 'gold', 'magenta', 'cyan', 'brown', 'orange')

# read in Casal2 MPD results
cas2_mpd <- list()
cas2_corr <- list()
cas2_conv <- list()
for (c in 1:num_C2_models)
{
    cas2_mpd_filename <- file.path(base_dir, 'Casal2', C2_subdir[c], mpd_run)
    cas2_mpd[[c]] <- casal2::extract.mpd(cas2_mpd_filename)

    cas2_corr[[c]] <- casal2::extract.correlation.matrix(mpd_out, file.path(base_dir, 'Casal2', C2_subdir[c]))

    cas2_conv[[c]] <- get_convergence_information(cas2_mpd_filename)
}

1.3 Tables

Tables of parameter estimates and objective function components for the CASAL and Casal2 model MPD results

CASAL parameter estimates
rownames Base_Model Sensitivity_1 Percent_Diff
q[CSacous].q                0.43115 0.43116 -0.003
q[WCacous].q                0.30025 0.30024 0.003
q[CRsum].q                0.07431 0.07431 0.007
q[SAsum].q                0.05925 0.05925 0.001
q[SAaut].q                0.06717 0.06717 0.002
relative_abundance[CRsumbio].cv_process_error 0.14048 0.14048 0.001
relative_abundance[SAsumbio].cv_process_error 0.43061 0.43061 -0.000
log_B0_total                  14.33950 14.33950 0.000
B0_prop_stock1                  0.32872 0.32871 0.003
natural_mortality.all              0.29598 0.29598 -0.002
selectivity[Wspsl].shift_a          -0.21743 -0.21744 -0.005
recruitment[E].YCS1             0.44781 0.44782 -0.001
recruitment[E].YCS2             0.72678 0.72680 -0.002
recruitment[E].YCS3             0.78452 0.78454 -0.002
recruitment[E].YCS4             1.35998 1.36001 -0.002
recruitment[E].YCS5             0.98774 0.98776 -0.002
recruitment[E].YCS6             0.32931 0.32932 -0.001
recruitment[E].YCS7             0.44722 0.44723 -0.001
recruitment[E].YCS8             0.26556 0.26556 -0.001
recruitment[E].YCS9             0.43418 0.43418 0.000
recruitment[E].YCS10            0.48804 0.48804 -0.001
recruitment[E].YCS11            0.39569 0.39570 -0.000
recruitment[E].YCS12            0.21791 0.21791 -0.000
recruitment[E].YCS13            1.03641 1.03640 0.001
recruitment[E].YCS14            0.95892 0.95891 0.001
recruitment[E].YCS15            0.42899 0.42899 -0.001
recruitment[E].YCS16            0.13810 0.13810 -0.001
recruitment[E].YCS17            0.37501 0.37500 0.002
recruitment[E].YCS18            0.39166 0.39164 0.005
recruitment[E].YCS19            0.38310 0.38310 0.000
recruitment[E].YCS20            0.56286 0.56285 0.002
recruitment[E].YCS21            0.26399 0.26399 0.000
recruitment[E].YCS22            0.28447 0.28447 0.000
recruitment[E].YCS23            0.45004 0.45003 0.001
recruitment[E].YCS24            0.25308 0.25307 0.001
recruitment[E].YCS25            0.06000 0.06000 0.000
recruitment[E].YCS26            0.46220 0.46220 0.001
recruitment[E].YCS27            0.16165 0.16165 -0.001
recruitment[E].YCS28            0.26690 0.26689 0.003
recruitment[E].YCS29            0.50038 0.50038 0.000
recruitment[E].YCS30            0.52003 0.52003 0.001
recruitment[E].YCS31            0.50883 0.50882 0.000
recruitment[E].YCS32            0.60056 0.60055 0.001
recruitment[E].YCS33            0.52215 0.52215 0.000
recruitment[E].YCS34            0.37302 0.37300 0.003
recruitment[E].YCS35            0.52681 0.52681 -0.000
recruitment[E].YCS36            0.15864 0.15864 0.001
recruitment[E].YCS37            0.66058 0.66058 0.000
recruitment[E].YCS38            0.15599 0.15600 -0.001
recruitment[E].YCS39            0.18584 0.18584 -0.001
recruitment[E].YCS40            0.73063 0.73062 0.001
recruitment[E].YCS41            0.73906 0.73900 0.008
recruitment[E].YCS42            1.19514 1.19515 -0.001
recruitment[W].YCS1             0.45702 0.45702 -0.000
recruitment[W].YCS2             0.82218 0.82218 0.000
recruitment[W].YCS3             0.67001 0.67000 0.000
recruitment[W].YCS4             0.65390 0.65390 0.000
recruitment[W].YCS5             0.46590 0.46590 0.001
recruitment[W].YCS6             0.51512 0.51512 -0.000
recruitment[W].YCS7             0.57824 0.57824 0.000
recruitment[W].YCS8             0.47328 0.47328 0.001
recruitment[W].YCS9             0.87191 0.87192 -0.002
recruitment[W].YCS10            0.83644 0.83646 -0.003
recruitment[W].YCS11            0.32024 0.32024 -0.002
recruitment[W].YCS12            0.16497 0.16497 0.001
recruitment[W].YCS13            1.47554 1.47555 -0.001
recruitment[W].YCS14            0.77274 0.77276 -0.002
recruitment[W].YCS15            0.08135 0.08134 0.000
recruitment[W].YCS16            0.08158 0.08158 0.001
recruitment[W].YCS17            0.89280 0.89280 0.000
recruitment[W].YCS18            1.26777 1.26778 -0.001
recruitment[W].YCS19            0.37722 0.37722 0.000
recruitment[W].YCS20            0.95940 0.95939 0.000
recruitment[W].YCS21            0.22113 0.22113 0.001
recruitment[W].YCS22            0.19718 0.19718 0.001
recruitment[W].YCS23            0.23794 0.23795 -0.003
recruitment[W].YCS24            0.20291 0.20292 -0.002
recruitment[W].YCS25            0.06000 0.06000 0.000
recruitment[W].YCS26            0.37304 0.37304 -0.002
recruitment[W].YCS27            0.09548 0.09548 -0.001
recruitment[W].YCS28            0.79558 0.79558 -0.000
recruitment[W].YCS29            0.46516 0.46517 -0.001
recruitment[W].YCS30            0.51662 0.51662 -0.001
recruitment[W].YCS31            0.51126 0.51126 -0.000
recruitment[W].YCS32            0.75544 0.75544 0.000
recruitment[W].YCS33            0.63892 0.63891 0.001
recruitment[W].YCS34            0.61686 0.61686 -0.001
recruitment[W].YCS35            0.57841 0.57841 0.001
recruitment[W].YCS36            0.14547 0.14547 -0.001
recruitment[W].YCS37            0.97479 0.97478 0.001
recruitment[W].YCS38            0.16016 0.16016 0.002
recruitment[W].YCS39            0.24444 0.24444 0.001
recruitment[W].YCS40            2.19661 2.19657 0.002
recruitment[W].YCS41            2.15201 2.15204 -0.001
recruitment[W].YCS42            1.19538 1.19539 -0.001
migration[Whome].rates_all1       0.32784 0.32783 0.003
migration[Whome].rates_all2       0.30780 0.30778 0.006
migration[Whome].rates_all3       0.47624 0.47620 0.007
migration[Whome].rates_all4       0.14184 0.14184 0.001
migration[Whome].rates_all5       0.46851 0.46836 0.033
migration[Whome].rates_all6       0.74978 0.75026 -0.064
migration[Whome].rates_all7       0.79469 0.78841 0.790
migration[Whome].rates_all8       1.00000 1.00000 0.000
migration[Espmg].rates_all1       0.00566 0.00566 0.002
migration[Espmg].rates_all2       0.16021 0.16021 0.002
migration[Espmg].rates_all3       0.32412 0.32411 0.002
migration[Espmg].rates_all4       0.43124 0.43123 0.002
migration[Espmg].rates_all5       0.38546 0.38546 0.002
migration[Espmg].rates_all6       0.46536 0.46535 0.002
migration[Espmg].rates_all7       0.59902 0.59902 0.001
migration[Espmg].rates_all8       0.60778 0.60778 0.001
migration[Wspmg].rates_all1       0.10614 0.10614 0.000
migration[Wspmg].rates_all2       0.28774 0.28774 -0.001
migration[Wspmg].rates_all3       0.27076 0.27076 0.000
migration[Wspmg].rates_all4       0.45525 0.45525 0.000
migration[Wspmg].rates_all5       0.48465 0.48466 -0.001
migration[Wspmg].rates_all6       0.52713 0.52711 0.005
migration[Wspmg].rates_all7       0.66056 0.66062 -0.009
migration[Wspmg].rates_all8       0.59628 0.59627 0.001
selectivity[Enspsl].all1             3.56717 3.56716 0.000
selectivity[Enspsl].all2             1.22109 1.22108 0.001
selectivity[Enspsl].all3             15.71130 15.71370 -0.015
selectivity[Wnspsl].all1             5.34113 5.34117 -0.001
selectivity[Wnspsl].all2             2.83932 2.83934 -0.001
selectivity[Wnspsl].all3             13.37880 13.37890 -0.001
selectivity[Espsl].all1             10.87880 10.87880 0.000
selectivity[Espsl].all2             11.67780 11.67750 0.003
selectivity[Espsl].all3             3.98159 3.98165 -0.002
selectivity[CRsl].all1             1.00000 1.00000 0.000
selectivity[CRsl].all2             1.00000 1.00000 0.000
selectivity[CRsl].all3             9.30026 9.30042 -0.002
selectivity[SAsl].all1             14.25200 14.25230 -0.002
selectivity[SAsl].all2             8.75543 8.75558 -0.002
selectivity[SAsl].all3             1.31548 1.31519 0.022
Casal2 parameter estimates
rownames betadiff_casal_flags_on betadiff_casal_flags_off betadiff_casal_flags_on_low_tol cppad_casal_flags_on cppad_casal_flags_off
catchability[CSacous].q                0.428 0.427 0.431 0.431 0.431
catchability[WCacous].q                0.312 0.313 0.300 0.300 0.300
catchability[CRsum].q                0.068 0.068 0.074 0.074 0.074
catchability[SAsum].q                0.067 0.067 0.059 0.059 0.059
catchability[SAaut].q                0.076 0.076 0.067 0.067 0.067
observation[CRsumbio].process_error    0.141 0.141 0.140 0.140 0.140
observation[SAsumbio].process_error    0.424 0.424 0.431 0.431 0.431
process[recruit_E].b0               547497.000 548395.000 555125.000 555125.000 555125.000
process[recruit_W].b0               1083360.000 1082670.000 1133680.000 1133680.000 1133680.000
process[Instant_mortality].m{west.sa}              0.313 0.313 0.296 0.296 0.296
time_varying[shifted_mu].a                -0.230 -0.229 -0.217 -0.217 -0.217
process[recruit_E].ycs_values{1975} 0.446 0.445 0.448 0.448 0.448
process[recruit_E].ycs_values{1976} 0.731 0.729 0.727 0.727 0.727
process[recruit_E].ycs_values{1977} 0.793 0.792 0.785 0.785 0.785
process[recruit_E].ycs_values{1978} 1.388 1.386 1.360 1.360 1.360
process[recruit_E].ycs_values{1979} 1.004 1.003 0.988 0.988 0.988
process[recruit_E].ycs_values{1980} 0.335 0.334 0.329 0.329 0.329
process[recruit_E].ycs_values{1981} 0.455 0.455 0.447 0.447 0.447
process[recruit_E].ycs_values{1982} 0.270 0.270 0.266 0.266 0.266
process[recruit_E].ycs_values{1983} 0.441 0.441 0.434 0.434 0.434
process[recruit_E].ycs_values{1984} 0.494 0.494 0.488 0.488 0.488
process[recruit_E].ycs_values{1985} 0.400 0.400 0.396 0.396 0.396
process[recruit_E].ycs_values{1986} 0.220 0.220 0.218 0.218 0.218
process[recruit_E].ycs_values{1987} 1.041 1.041 1.036 1.036 1.036
process[recruit_E].ycs_values{1988} 0.963 0.963 0.959 0.959 0.959
process[recruit_E].ycs_values{1989} 0.431 0.431 0.429 0.429 0.429
process[recruit_E].ycs_values{1990} 0.138 0.138 0.138 0.138 0.138
process[recruit_E].ycs_values{1991} 0.373 0.374 0.375 0.375 0.375
process[recruit_E].ycs_values{1992} 0.388 0.389 0.392 0.392 0.392
process[recruit_E].ycs_values{1993} 0.380 0.380 0.383 0.383 0.383
process[recruit_E].ycs_values{1994} 0.557 0.558 0.563 0.563 0.563
process[recruit_E].ycs_values{1995} 0.262 0.262 0.264 0.264 0.264
process[recruit_E].ycs_values{1996} 0.282 0.283 0.284 0.284 0.284
process[recruit_E].ycs_values{1997} 0.448 0.448 0.450 0.450 0.450
process[recruit_E].ycs_values{1998} 0.252 0.253 0.253 0.253 0.253
process[recruit_E].ycs_values{1999} 0.060 0.060 0.060 0.060 0.060
process[recruit_E].ycs_values{2000} 0.460 0.460 0.462 0.462 0.462
process[recruit_E].ycs_values{2001} 0.161 0.161 0.162 0.162 0.162
process[recruit_E].ycs_values{2002} 0.266 0.266 0.267 0.267 0.267
process[recruit_E].ycs_values{2003} 0.499 0.500 0.500 0.500 0.500
process[recruit_E].ycs_values{2004} 0.520 0.521 0.520 0.520 0.520
process[recruit_E].ycs_values{2005} 0.508 0.508 0.509 0.509 0.509
process[recruit_E].ycs_values{2006} 0.600 0.600 0.601 0.601 0.601
process[recruit_E].ycs_values{2007} 0.521 0.521 0.522 0.522 0.522
process[recruit_E].ycs_values{2008} 0.370 0.371 0.373 0.373 0.373
process[recruit_E].ycs_values{2009} 0.528 0.528 0.527 0.527 0.527
process[recruit_E].ycs_values{2010} 0.157 0.157 0.159 0.159 0.159
process[recruit_E].ycs_values{2011} 0.660 0.660 0.661 0.661 0.661
process[recruit_E].ycs_values{2012} 0.156 0.156 0.156 0.156 0.156
process[recruit_E].ycs_values{2013} 0.185 0.185 0.186 0.186 0.186
process[recruit_E].ycs_values{2014} 0.728 0.728 0.731 0.731 0.731
process[recruit_E].ycs_values{2015} 0.752 0.756 0.739 0.739 0.739
process[recruit_E].ycs_values{2016} 1.189 1.189 1.195 1.195 1.195
process[recruit_W].ycs_values{1975} 0.459 0.459 0.457 0.457 0.457
process[recruit_W].ycs_values{1976} 0.814 0.814 0.822 0.822 0.822
process[recruit_W].ycs_values{1977} 0.685 0.685 0.670 0.670 0.670
process[recruit_W].ycs_values{1978} 0.669 0.669 0.654 0.654 0.654
process[recruit_W].ycs_values{1979} 0.473 0.473 0.466 0.466 0.466
process[recruit_W].ycs_values{1980} 0.521 0.521 0.515 0.515 0.515
process[recruit_W].ycs_values{1981} 0.583 0.583 0.578 0.578 0.578
process[recruit_W].ycs_values{1982} 0.474 0.475 0.473 0.473 0.473
process[recruit_W].ycs_values{1983} 0.874 0.873 0.872 0.872 0.872
process[recruit_W].ycs_values{1984} 0.837 0.836 0.836 0.836 0.836
process[recruit_W].ycs_values{1985} 0.320 0.319 0.320 0.320 0.320
process[recruit_W].ycs_values{1986} 0.165 0.165 0.165 0.165 0.165
process[recruit_W].ycs_values{1987} 1.481 1.480 1.476 1.476 1.476
process[recruit_W].ycs_values{1988} 0.775 0.775 0.773 0.773 0.773
process[recruit_W].ycs_values{1989} 0.081 0.081 0.081 0.081 0.081
process[recruit_W].ycs_values{1990} 0.082 0.082 0.082 0.082 0.082
process[recruit_W].ycs_values{1991} 0.892 0.892 0.893 0.893 0.893
process[recruit_W].ycs_values{1992} 1.268 1.268 1.268 1.268 1.268
process[recruit_W].ycs_values{1993} 0.376 0.377 0.377 0.377 0.377
process[recruit_W].ycs_values{1994} 0.956 0.956 0.959 0.959 0.959
process[recruit_W].ycs_values{1995} 0.220 0.220 0.221 0.221 0.221
process[recruit_W].ycs_values{1996} 0.196 0.196 0.197 0.197 0.197
process[recruit_W].ycs_values{1997} 0.237 0.236 0.238 0.238 0.238
process[recruit_W].ycs_values{1998} 0.201 0.200 0.203 0.203 0.203
process[recruit_W].ycs_values{1999} 0.060 0.060 0.060 0.060 0.060
process[recruit_W].ycs_values{2000} 0.371 0.371 0.373 0.373 0.373
process[recruit_W].ycs_values{2001} 0.095 0.095 0.095 0.095 0.095
process[recruit_W].ycs_values{2002} 0.796 0.796 0.796 0.796 0.796
process[recruit_W].ycs_values{2003} 0.466 0.466 0.465 0.465 0.465
process[recruit_W].ycs_values{2004} 0.518 0.518 0.517 0.517 0.517
process[recruit_W].ycs_values{2005} 0.513 0.513 0.511 0.511 0.511
process[recruit_W].ycs_values{2006} 0.755 0.755 0.755 0.755 0.755
process[recruit_W].ycs_values{2007} 0.637 0.637 0.639 0.639 0.639
process[recruit_W].ycs_values{2008} 0.616 0.616 0.617 0.617 0.617
process[recruit_W].ycs_values{2009} 0.573 0.573 0.578 0.578 0.578
process[recruit_W].ycs_values{2010} 0.146 0.146 0.145 0.145 0.145
process[recruit_W].ycs_values{2011} 0.967 0.967 0.975 0.975 0.975
process[recruit_W].ycs_values{2012} 0.159 0.159 0.160 0.160 0.160
process[recruit_W].ycs_values{2013} 0.244 0.244 0.244 0.244 0.244
process[recruit_W].ycs_values{2014} 2.185 2.187 2.197 2.197 2.197
process[recruit_W].ycs_values{2015} 2.150 2.148 2.152 2.152 2.152
process[recruit_W].ycs_values{2016} 1.190 1.189 1.195 1.195 1.195
selectivity[sel_Whome].v{1}             0.321 0.322 0.328 0.328 0.328
selectivity[sel_Whome].v{2}             0.284 0.285 0.308 0.308 0.308
selectivity[sel_Whome].v{3}             0.458 0.460 0.476 0.476 0.476
selectivity[sel_Whome].v{4}             0.135 0.135 0.142 0.142 0.142
selectivity[sel_Whome].v{5}             0.432 0.438 0.468 0.468 0.468
selectivity[sel_Whome].v{6}             0.719 0.708 0.750 0.750 0.750
selectivity[sel_Whome].v{7}             0.760 1.000 0.788 0.788 0.788
selectivity[sel_Espmg].v{1}             0.006 0.006 0.006 0.006 0.006
selectivity[sel_Espmg].v{2}             0.159 0.159 0.160 0.160 0.160
selectivity[sel_Espmg].v{3}             0.321 0.321 0.324 0.324 0.324
selectivity[sel_Espmg].v{4}             0.425 0.425 0.431 0.431 0.431
selectivity[sel_Espmg].v{5}             0.379 0.380 0.385 0.385 0.385
selectivity[sel_Espmg].v{6}             0.458 0.459 0.465 0.465 0.465
selectivity[sel_Espmg].v{7}             0.590 0.591 0.599 0.599 0.599
selectivity[sel_Espmg].v{8}             0.611 0.611 0.608 0.608 0.608
selectivity[sel_Wspmg].v{1}             0.106 0.106 0.106 0.106 0.106
selectivity[sel_Wspmg].v{2}             0.291 0.291 0.288 0.288 0.288
selectivity[sel_Wspmg].v{3}             0.271 0.271 0.271 0.271 0.271
selectivity[sel_Wspmg].v{4}             0.455 0.455 0.455 0.455 0.455
selectivity[sel_Wspmg].v{5}             0.484 0.484 0.485 0.485 0.485
selectivity[sel_Wspmg].v{6}             0.524 0.525 0.527 0.527 0.527
selectivity[sel_Wspmg].v{7}             0.656 0.653 0.661 0.661 0.661
selectivity[sel_Wspmg].v{8}             0.601 0.602 0.596 0.596 0.596
selectivity[Enspsl].mu               3.572 3.573 3.567 3.567 3.567
selectivity[Enspsl].sigma_l          1.222 1.222 1.221 1.221 1.221
selectivity[Enspsl].sigma_r          50.000 50.000 15.714 15.714 15.714
selectivity[Wnspsl].mu               5.330 5.328 5.341 5.341 5.341
selectivity[Wnspsl].sigma_l          2.778 2.778 2.839 2.839 2.839
selectivity[Wnspsl].sigma_r          49.998 50.000 13.379 13.379 13.379
selectivity[Espsl].mu               11.255 11.257 10.879 10.879 10.879
selectivity[Espsl].sigma_l          10.859 10.867 11.678 11.678 11.678
selectivity[Espsl].sigma_r          4.029 4.026 3.982 3.982 3.982
selectivity[CRsl].mu               1.000 1.000 1.000 1.000 1.000
selectivity[CRsl].sigma_l          46.483 46.482 45.986 25.500 25.500
selectivity[CRsl].sigma_r          10.088 10.093 9.300 9.300 9.300
selectivity[SAsl].mu               14.865 14.868 14.252 14.252 14.252
selectivity[SAsl].sigma_l          8.693 8.695 8.756 8.756 8.756
selectivity[SAsl].sigma_r          1.000 1.000 1.315 1.315 1.315
Casal2 parameter estimates: Percent Difference from betadiff_casal_flags_on
rownames betadiff_casal_flags_on betadiff_casal_flags_off betadiff_casal_flags_on_low_tol cppad_casal_flags_on cppad_casal_flags_off
catchability[CSacous].q                0.000 0.145 -0.772 -0.772 -0.772
catchability[WCacous].q                0.000 -0.128 3.904 3.904 3.904
catchability[CRsum].q                0.000 -0.299 -9.090 -9.090 -9.090
catchability[SAsum].q                0.000 -0.172 11.645 11.645 11.645
catchability[SAaut].q                0.000 -0.189 11.196 11.196 11.196
observation[CRsumbio].process_error    0.000 -0.035 0.292 0.292 0.292
observation[SAsumbio].process_error    0.000 0.041 -1.604 -1.604 -1.604
process[recruit_E].b0               0.000 -0.164 -1.393 -1.393 -1.393
process[recruit_W].b0               0.000 0.064 -4.645 -4.645 -4.645
process[Instant_mortality].m{west.sa}              0.000 0.070 5.415 5.415 5.415
time_varying[shifted_mu].a                0.000 0.166 5.293 5.293 5.293
process[recruit_E].ycs_values{1975} 0.000 0.117 -0.431 -0.431 -0.431
process[recruit_E].ycs_values{1976} 0.000 0.156 0.510 0.510 0.510
process[recruit_E].ycs_values{1977} 0.000 0.127 1.107 1.107 1.107
process[recruit_E].ycs_values{1978} 0.000 0.131 2.036 2.036 2.036
process[recruit_E].ycs_values{1979} 0.000 0.107 1.590 1.590 1.590
process[recruit_E].ycs_values{1980} 0.000 0.073 1.584 1.584 1.584
process[recruit_E].ycs_values{1981} 0.000 0.072 1.783 1.783 1.783
process[recruit_E].ycs_values{1982} 0.000 0.055 1.605 1.605 1.605
process[recruit_E].ycs_values{1983} 0.000 -0.013 1.547 1.547 1.547
process[recruit_E].ycs_values{1984} 0.000 0.022 1.251 1.251 1.251
process[recruit_E].ycs_values{1985} 0.000 0.010 1.111 1.111 1.111
process[recruit_E].ycs_values{1986} 0.000 0.022 0.872 0.872 0.872
process[recruit_E].ycs_values{1987} 0.000 -0.035 0.437 0.437 0.437
process[recruit_E].ycs_values{1988} 0.000 -0.032 0.440 0.440 0.440
process[recruit_E].ycs_values{1989} 0.000 0.040 0.495 0.495 0.495
process[recruit_E].ycs_values{1990} 0.000 0.009 -0.091 -0.091 -0.091
process[recruit_E].ycs_values{1991} 0.000 -0.136 -0.443 -0.443 -0.443
process[recruit_E].ycs_values{1992} 0.000 -0.260 -1.050 -1.050 -1.050
process[recruit_E].ycs_values{1993} 0.000 -0.017 -0.847 -0.847 -0.847
process[recruit_E].ycs_values{1994} 0.000 -0.122 -1.003 -1.003 -1.003
process[recruit_E].ycs_values{1995} 0.000 0.001 -0.842 -0.842 -0.842
process[recruit_E].ycs_values{1996} 0.000 -0.008 -0.704 -0.704 -0.704
process[recruit_E].ycs_values{1997} 0.000 -0.049 -0.545 -0.545 -0.545
process[recruit_E].ycs_values{1998} 0.000 -0.061 -0.232 -0.232 -0.232
process[recruit_E].ycs_values{1999} 0.000 0.000 0.000 0.000 0.000
process[recruit_E].ycs_values{2000} 0.000 -0.057 -0.506 -0.506 -0.506
process[recruit_E].ycs_values{2001} 0.000 0.007 -0.482 -0.482 -0.482
process[recruit_E].ycs_values{2002} 0.000 -0.191 -0.486 -0.486 -0.486
process[recruit_E].ycs_values{2003} 0.000 -0.031 -0.190 -0.190 -0.190
process[recruit_E].ycs_values{2004} 0.000 -0.043 0.054 0.054 0.054
process[recruit_E].ycs_values{2005} 0.000 -0.021 -0.254 -0.254 -0.254
process[recruit_E].ycs_values{2006} 0.000 -0.046 -0.116 -0.116 -0.116
process[recruit_E].ycs_values{2007} 0.000 -0.015 -0.259 -0.259 -0.259
process[recruit_E].ycs_values{2008} 0.000 -0.184 -0.685 -0.685 -0.685
process[recruit_E].ycs_values{2009} 0.000 -0.007 0.199 0.199 0.199
process[recruit_E].ycs_values{2010} 0.000 -0.033 -1.107 -1.107 -1.107
process[recruit_E].ycs_values{2011} 0.000 -0.027 -0.044 -0.044 -0.044
process[recruit_E].ycs_values{2012} 0.000 0.030 -0.264 -0.264 -0.264
process[recruit_E].ycs_values{2013} 0.000 -0.001 -0.357 -0.357 -0.357
process[recruit_E].ycs_values{2014} 0.000 -0.049 -0.349 -0.349 -0.349
process[recruit_E].ycs_values{2015} 0.000 -0.472 1.769 1.769 1.769
process[recruit_E].ycs_values{2016} 0.000 0.045 -0.484 -0.484 -0.484
process[recruit_W].ycs_values{1975} 0.000 0.005 0.371 0.371 0.371
process[recruit_W].ycs_values{1976} 0.000 0.014 -0.963 -0.963 -0.963
process[recruit_W].ycs_values{1977} 0.000 -0.030 2.180 2.180 2.180
process[recruit_W].ycs_values{1978} 0.000 -0.047 2.265 2.265 2.265
process[recruit_W].ycs_values{1979} 0.000 -0.067 1.451 1.451 1.451
process[recruit_W].ycs_values{1980} 0.000 -0.023 1.038 1.038 1.038
process[recruit_W].ycs_values{1981} 0.000 -0.027 0.865 0.865 0.865
process[recruit_W].ycs_values{1982} 0.000 -0.047 0.237 0.237 0.237
process[recruit_W].ycs_values{1983} 0.000 0.076 0.209 0.209 0.209
process[recruit_W].ycs_values{1984} 0.000 0.125 0.109 0.109 0.109
process[recruit_W].ycs_values{1985} 0.000 0.103 -0.153 -0.153 -0.153
process[recruit_W].ycs_values{1986} 0.000 -0.055 -0.076 -0.076 -0.076
process[recruit_W].ycs_values{1987} 0.000 0.040 0.345 0.345 0.345
process[recruit_W].ycs_values{1988} 0.000 0.098 0.327 0.327 0.327
process[recruit_W].ycs_values{1989} 0.000 -0.003 0.158 0.158 0.158
process[recruit_W].ycs_values{1990} 0.000 -0.037 0.267 0.267 0.267
process[recruit_W].ycs_values{1991} 0.000 -0.010 -0.102 -0.102 -0.102
process[recruit_W].ycs_values{1992} 0.000 0.021 0.004 0.004 0.004
process[recruit_W].ycs_values{1993} 0.000 -0.011 -0.196 -0.196 -0.196
process[recruit_W].ycs_values{1994} 0.000 -0.024 -0.400 -0.400 -0.400
process[recruit_W].ycs_values{1995} 0.000 -0.053 -0.446 -0.446 -0.446
process[recruit_W].ycs_values{1996} 0.000 -0.029 -0.468 -0.468 -0.468
process[recruit_W].ycs_values{1997} 0.000 0.147 -0.496 -0.496 -0.496
process[recruit_W].ycs_values{1998} 0.000 0.099 -1.151 -1.151 -1.151
process[recruit_W].ycs_values{1999} 0.000 0.000 0.000 0.000 0.000
process[recruit_W].ycs_values{2000} 0.000 0.117 -0.463 -0.463 -0.463
process[recruit_W].ycs_values{2001} 0.000 0.037 -0.409 -0.409 -0.409
process[recruit_W].ycs_values{2002} 0.000 0.029 0.022 0.022 0.022
process[recruit_W].ycs_values{2003} 0.000 0.038 0.153 0.153 0.153
process[recruit_W].ycs_values{2004} 0.000 0.045 0.235 0.235 0.235
process[recruit_W].ycs_values{2005} 0.000 0.016 0.366 0.366 0.366
process[recruit_W].ycs_values{2006} 0.000 0.015 -0.076 -0.076 -0.076
process[recruit_W].ycs_values{2007} 0.000 -0.025 -0.338 -0.338 -0.338
process[recruit_W].ycs_values{2008} 0.000 0.080 -0.118 -0.118 -0.118
process[recruit_W].ycs_values{2009} 0.000 -0.024 -0.958 -0.958 -0.958
process[recruit_W].ycs_values{2010} 0.000 0.055 0.538 0.538 0.538
process[recruit_W].ycs_values{2011} 0.000 -0.041 -0.837 -0.837 -0.837
process[recruit_W].ycs_values{2012} 0.000 -0.075 -0.750 -0.750 -0.750
process[recruit_W].ycs_values{2013} 0.000 -0.052 -0.224 -0.224 -0.224
process[recruit_W].ycs_values{2014} 0.000 -0.081 -0.541 -0.541 -0.541
process[recruit_W].ycs_values{2015} 0.000 0.098 -0.109 -0.109 -0.109
process[recruit_W].ycs_values{2016} 0.000 0.045 -0.485 -0.485 -0.485
selectivity[sel_Whome].v{1}             0.000 -0.177 -2.014 -2.014 -2.014
selectivity[sel_Whome].v{2}             0.000 -0.372 -8.428 -8.428 -8.428
selectivity[sel_Whome].v{3}             0.000 -0.431 -4.031 -4.031 -4.031
selectivity[sel_Whome].v{4}             0.000 -0.344 -5.047 -5.047 -5.047
selectivity[sel_Whome].v{5}             0.000 -1.419 -8.494 -8.494 -8.494
selectivity[sel_Whome].v{6}             0.000 1.487 -4.361 -4.361 -4.361
selectivity[sel_Whome].v{7}             0.000 -31.631 -3.792 -3.792 -3.792
selectivity[sel_Espmg].v{1}             0.000 -0.081 -0.412 -0.412 -0.412
selectivity[sel_Espmg].v{2}             0.000 -0.074 -0.757 -0.757 -0.757
selectivity[sel_Espmg].v{3}             0.000 -0.079 -1.107 -1.107 -1.107
selectivity[sel_Espmg].v{4}             0.000 -0.079 -1.434 -1.434 -1.434
selectivity[sel_Espmg].v{5}             0.000 -0.069 -1.633 -1.633 -1.633
selectivity[sel_Espmg].v{6}             0.000 -0.063 -1.500 -1.500 -1.500
selectivity[sel_Espmg].v{7}             0.000 -0.062 -1.452 -1.452 -1.452
selectivity[sel_Espmg].v{8}             0.000 -0.060 0.468 0.468 0.468
selectivity[sel_Wspmg].v{1}             0.000 -0.033 -0.586 -0.586 -0.586
selectivity[sel_Wspmg].v{2}             0.000 0.010 1.252 1.252 1.252
selectivity[sel_Wspmg].v{3}             0.000 -0.008 -0.037 -0.037 -0.037
selectivity[sel_Wspmg].v{4}             0.000 -0.016 -0.012 -0.012 -0.012
selectivity[sel_Wspmg].v{5}             0.000 0.026 -0.066 -0.066 -0.066
selectivity[sel_Wspmg].v{6}             0.000 -0.187 -0.597 -0.597 -0.597
selectivity[sel_Wspmg].v{7}             0.000 0.447 -0.718 -0.718 -0.718
selectivity[sel_Wspmg].v{8}             0.000 -0.088 0.788 0.788 0.788
selectivity[Enspsl].mu               0.000 -0.025 0.143 0.143 0.143
selectivity[Enspsl].sigma_l          0.000 -0.031 0.080 0.080 0.080
selectivity[Enspsl].sigma_r          0.000 -0.000 68.573 68.573 68.573
selectivity[Wnspsl].mu               0.000 0.021 -0.219 -0.219 -0.219
selectivity[Wnspsl].sigma_l          0.000 0.011 -2.191 -2.191 -2.191
selectivity[Wnspsl].sigma_r          0.000 -0.003 73.241 73.241 73.241
selectivity[Espsl].mu               0.000 -0.019 3.339 3.339 3.339
selectivity[Espsl].sigma_l          0.000 -0.073 -7.534 -7.534 -7.534
selectivity[Espsl].sigma_r          0.000 0.071 1.178 1.178 1.178
selectivity[CRsl].mu               0.000 0.000 0.000 0.000 0.000
selectivity[CRsl].sigma_l          0.000 0.002 1.069 45.141 45.141
selectivity[CRsl].sigma_r          0.000 -0.050 7.805 7.805 7.805
selectivity[SAsl].mu               0.000 -0.025 4.120 4.120 4.120
selectivity[SAsl].sigma_l          0.000 -0.022 -0.722 -0.722 -0.722
selectivity[SAsl].sigma_r          0.000 0.000 -31.519 -31.519 -31.519

CASAL objective function component values
Component Base_Model Sensitivity_1
CRsumbio -31.032 -31.032
CSacous -9.835 -9.836
SAautbio -3.984 -3.984
SAsumbio -6.146 -6.146
WCacous -5.916 -5.916
CRsumage 323.850 323.850
SAautage 35.117 35.118
SAsumage 221.149 221.150
EnspOLF 24.276 24.276
Enspage 256.069 256.068
Espage 581.502 581.502
WnspOLF 70.250 70.250
Wnspage 157.004 157.004
Wspage 364.751 364.751
pspawn -13.129 -13.129
pspawn_1993 -6.015 -6.015
prior_on_q[CSacous].q -0.839 -0.839
prior_on_q[WCacous].q -1.202 -1.202
prior_on_q[CRsum].q -2.207 -2.207
prior_on_q[SAsum].q -1.557 -1.557
prior_on_q[SAaut].q -1.762 -1.762
prior_on_relative_abundance[CRsumbio].cv_process_error 0.000 0.000
prior_on_relative_abundance[SAsumbio].cv_process_error 0.000 0.000
prior_on_log_B0_total 0.000 0.000
prior_on_B0_prop_stock1 12.555 12.555
prior_on_natural_mortality.all -1.217 -1.217
prior_on_selectivity[Wspsl].shift_a 0.378 0.378
prior_on_recruitment[E].YCS -13.750 -13.750
prior_on_recruitment[W].YCS -1.785 -1.785
prior_on_migration[Whome].rates_all 0.000 0.000
prior_on_migration[Espmg].rates_all 0.000 0.000
prior_on_migration[Wspmg].rates_all 0.000 0.000
prior_on_selectivity[Enspsl].all 0.000 0.000
prior_on_selectivity[Wnspsl].all 0.000 0.000
prior_on_selectivity[Espsl].all 0.000 0.000
prior_on_selectivity[CRsl].all 0.000 0.000
prior_on_selectivity[SAsl].all 0.000 0.000
Umax_Ensp1 0.000 0.000
Umax_Ensp2 0.000 0.000
Umax_Wnsp1 0.000 0.000
Umax_Wnsp2 0.000 0.000
Umax_Esp 0.000 0.000
Umax_Wsp 0.000 0.000
YCS.eq.16 0.000 0.000
sp.migr 0.013 0.013
Total 1946.540 1946.540
Casal2 objective function component values
rownames betadiff_casal_flags_on betadiff_casal_flags_off betadiff_casal_flags_on_low_tol cppad_casal_flags_on cppad_casal_flags_off
observation->CSacous-1991                  2.344 2.342 2.404 2.404 2.404
observation->CSacous-1993                  -0.577 -0.577 -0.588 -0.588 -0.588
observation->CSacous-1994                  -0.016 -0.016 -0.030 -0.030 -0.030
observation->CSacous-1995                  -0.562 -0.562 -0.567 -0.567 -0.567
observation->CSacous-1996                  -0.027 -0.027 0.013 0.013 0.013
observation->CSacous-1997                  -0.912 -0.912 -0.899 -0.899 -0.899
observation->CSacous-1998                  -0.055 -0.054 -0.004 -0.004 -0.004
observation->CSacous-1999                  -1.025 -1.025 -1.012 -1.012 -1.012
observation->CSacous-2001                  -1.217 -1.217 -1.211 -1.211 -1.211
observation->CSacous-2002                  -0.297 -0.297 -0.338 -0.338 -0.338
observation->CSacous-2003                  -0.845 -0.845 -0.863 -0.863 -0.863
observation->CSacous-2005                  -0.929 -0.929 -0.927 -0.927 -0.927
observation->CSacous-2006                  -0.939 -0.939 -0.950 -0.950 -0.950
observation->CSacous-2007                  -0.623 -0.623 -0.600 -0.600 -0.600
observation->CSacous-2008                  -1.094 -1.094 -1.125 -1.125 -1.125
observation->CSacous-2009                  -0.186 -0.185 -0.110 -0.110 -0.110
observation->CSacous-2011                  -1.021 -1.021 -1.000 -1.000 -1.000
observation->CSacous-2013                  -0.892 -0.892 -0.855 -0.855 -0.855
observation->CSacous-2015                  -0.191 -0.190 -0.150 -0.150 -0.150
observation->CSacous-2017                  -1.020 -1.020 -1.021 -1.021 -1.021
observation->WCacous-1988                  -0.588 -0.588 -0.588 -0.588 -0.588
observation->WCacous-1989                  -0.184 -0.183 -0.168 -0.168 -0.168
observation->WCacous-1990                  -0.331 -0.332 -0.327 -0.327 -0.327
observation->WCacous-1991                  -0.425 -0.425 -0.425 -0.425 -0.425
observation->WCacous-1992                  -0.744 -0.744 -0.744 -0.744 -0.744
observation->WCacous-1993                  -0.371 -0.369 -0.374 -0.374 -0.374
observation->WCacous-1997                  -0.055 -0.054 -0.043 -0.043 -0.043
observation->WCacous-2000                  -1.264 -1.265 -1.272 -1.272 -1.272
observation->WCacous-2012                  -0.906 -0.906 -0.896 -0.896 -0.896
observation->WCacous-2013                  -1.079 -1.079 -1.079 -1.079 -1.079
observation->CRsumbio-1992                  -1.634 -1.634 -1.647 -1.647 -1.647
observation->CRsumbio-1993                  0.041 0.040 0.043 0.043 0.043
observation->CRsumbio-1994                  -1.761 -1.760 -1.760 -1.760 -1.760
observation->CRsumbio-1995                  -1.431 -1.425 -1.406 -1.406 -1.406
observation->CRsumbio-1996                  -1.691 -1.690 -1.699 -1.699 -1.699
observation->CRsumbio-1997                  -0.609 -0.607 -0.635 -0.635 -0.635
observation->CRsumbio-1998                  -1.205 -1.203 -1.210 -1.210 -1.210
observation->CRsumbio-1999                  -0.950 -0.953 -0.919 -0.919 -0.919
observation->CRsumbio-2000                  -1.629 -1.629 -1.632 -1.632 -1.632
observation->CRsumbio-2001                  -1.773 -1.773 -1.775 -1.775 -1.775
observation->CRsumbio-2002                  -1.205 -1.209 -1.206 -1.206 -1.206
observation->CRsumbio-2003                  -1.748 -1.749 -1.744 -1.744 -1.744
observation->CRsumbio-2004                  -0.821 -0.818 -0.803 -0.803 -0.803
observation->CRsumbio-2005                  -1.470 -1.471 -1.477 -1.477 -1.477
observation->CRsumbio-2006                  -1.314 -1.314 -1.306 -1.306 -1.306
observation->CRsumbio-2007                  -0.726 -0.727 -0.756 -0.756 -0.756
observation->CRsumbio-2008                  -0.376 -0.373 -0.402 -0.402 -0.402
observation->CRsumbio-2009                  -0.587 -0.586 -0.559 -0.559 -0.559
observation->CRsumbio-2010                  -1.350 -1.349 -1.359 -1.359 -1.359
observation->CRsumbio-2011                  -1.094 -1.096 -1.105 -1.105 -1.105
observation->CRsumbio-2012                  -1.438 -1.439 -1.448 -1.448 -1.448
observation->CRsumbio-2013                  -1.494 -1.494 -1.491 -1.491 -1.491
observation->CRsumbio-2014                  -1.768 -1.768 -1.770 -1.770 -1.770
observation->CRsumbio-2016                  -1.411 -1.410 -1.411 -1.411 -1.411
observation->CRsumbio-2018                  -1.552 -1.552 -1.554 -1.554 -1.554
observation->SAsumbio-1992                  -0.596 -0.597 -0.588 -0.588 -0.588
observation->SAsumbio-1993                  -0.314 -0.314 -0.309 -0.309 -0.309
observation->SAsumbio-1994                  0.194 0.196 0.186 0.186 0.186
observation->SAsumbio-2001                  -0.834 -0.835 -0.822 -0.822 -0.822
observation->SAsumbio-2002                  -0.794 -0.794 -0.777 -0.777 -0.777
observation->SAsumbio-2003                  -0.836 -0.837 -0.825 -0.825 -0.825
observation->SAsumbio-2004                  0.735 0.732 0.734 0.734 0.734
observation->SAsumbio-2005                  -0.255 -0.258 -0.218 -0.218 -0.218
observation->SAsumbio-2006                  -0.598 -0.600 -0.558 -0.558 -0.558
observation->SAsumbio-2007                  0.637 0.636 0.654 0.654 0.654
observation->SAsumbio-2008                  -0.463 -0.461 -0.462 -0.462 -0.462
observation->SAsumbio-2009                  -0.584 -0.585 -0.583 -0.583 -0.583
observation->SAsumbio-2010                  -0.146 -0.143 -0.147 -0.147 -0.147
observation->SAsumbio-2012                  -0.848 -0.848 -0.835 -0.835 -0.835
observation->SAsumbio-2013                  -0.796 -0.797 -0.786 -0.786 -0.786
observation->SAsumbio-2015                  -0.201 -0.200 -0.189 -0.189 -0.189
observation->SAsumbio-2017                  -0.638 -0.638 -0.620 -0.620 -0.620
observation->SAautbio-1992                  -1.460 -1.459 -1.450 -1.450 -1.450
observation->SAautbio-1996                  -1.232 -1.233 -1.250 -1.250 -1.250
observation->SAautbio-1998                  -1.306 -1.308 -1.284 -1.284 -1.284
observation->Espage-1988                  19.790 19.793 19.640 19.640 19.640
observation->Espage-1989                  15.857 15.858 15.862 15.862 15.862
observation->Espage-1990                  21.486 21.488 21.446 21.446 21.446
observation->Espage-1991                  20.640 20.641 20.612 20.612 20.612
observation->Espage-1992                  24.092 24.089 24.020 24.020 24.020
observation->Espage-1993                  19.304 19.305 19.328 19.328 19.328
observation->Espage-1994                  20.522 20.528 20.503 20.503 20.503
observation->Espage-1995                  22.530 22.541 22.559 22.559 22.559
observation->Espage-1996                  18.590 18.590 18.577 18.577 18.577
observation->Espage-1997                  20.978 20.974 20.964 20.964 20.964
observation->Espage-1998                  20.579 20.587 20.611 20.611 20.611
observation->Espage-1999                  22.875 22.886 22.862 22.862 22.862
observation->Espage-2000                  23.321 23.321 23.314 23.314 23.314
observation->Espage-2001                  20.735 20.735 20.773 20.773 20.773
observation->Espage-2002                  26.243 26.248 26.287 26.287 26.287
observation->Espage-2003                  19.118 19.120 19.120 19.120 19.120
observation->Espage-2004                  22.698 22.692 22.642 22.642 22.642
observation->Espage-2005                  26.120 26.122 26.070 26.070 26.070
observation->Espage-2006                  23.245 23.242 23.149 23.149 23.149
observation->Espage-2007                  25.277 25.274 25.210 25.210 25.210
observation->Espage-2008                  22.166 22.168 22.129 22.129 22.129
observation->Espage-2009                  18.691 18.691 18.690 18.690 18.690
observation->Espage-2010                  21.856 21.855 21.851 21.851 21.851
observation->Espage-2014                  20.614 20.616 20.648 20.648 20.648
observation->Espage-2015                  21.820 21.823 21.870 21.870 21.870
observation->Espage-2016                  21.203 21.204 21.204 21.204 21.204
observation->Espage-2017                  21.578 21.575 21.563 21.563 21.563
observation->Wspage-1988                  15.729 15.724 15.753 15.753 15.753
observation->Wspage-1989                  11.558 11.555 11.566 11.566 11.566
observation->Wspage-1990                  11.750 11.752 11.790 11.790 11.790
observation->Wspage-1991                  13.037 13.036 13.051 13.051 13.051
observation->Wspage-1992                  10.179 10.178 10.170 10.170 10.170
observation->Wspage-1993                  10.571 10.569 10.579 10.579 10.579
observation->Wspage-1994                  9.918 9.918 9.924 9.924 9.924
observation->Wspage-1995                  13.399 13.403 13.396 13.396 13.396
observation->Wspage-1996                  11.841 11.840 11.844 11.844 11.844
observation->Wspage-1997                  11.233 11.230 11.243 11.243 11.243
observation->Wspage-1998                  11.761 11.761 11.762 11.762 11.762
observation->Wspage-1999                  11.976 11.975 11.981 11.981 11.981
observation->Wspage-2000                  11.256 11.256 11.260 11.260 11.260
observation->Wspage-2001                  12.181 12.179 12.201 12.201 12.201
observation->Wspage-2002                  9.732 9.731 9.745 9.745 9.745
observation->Wspage-2003                  9.555 9.555 9.556 9.556 9.556
observation->Wspage-2004                  14.342 14.340 14.304 14.304 14.304
observation->Wspage-2005                  14.934 14.932 14.977 14.977 14.977
observation->Wspage-2006                  13.817 13.817 13.872 13.872 13.872
observation->Wspage-2007                  11.500 11.498 11.516 11.516 11.516
observation->Wspage-2008                  12.196 12.195 12.194 12.194 12.194
observation->Wspage-2009                  12.635 12.631 12.594 12.594 12.594
observation->Wspage-2010                  11.317 11.316 11.299 11.299 11.299
observation->Wspage-2011                  11.155 11.154 11.137 11.137 11.137
observation->Wspage-2012                  11.866 11.867 11.860 11.860 11.860
observation->Wspage-2013                  11.043 11.041 11.065 11.065 11.065
observation->Wspage-2014                  12.944 12.944 12.941 12.941 12.941
observation->Wspage-2015                  13.136 13.134 13.101 13.101 13.101
observation->Wspage-2016                  14.268 14.272 14.273 14.273 14.273
observation->Wspage-2017                  13.746 13.745 13.797 13.797 13.797
observation->EnspOLF-1992                  4.396 4.397 4.388 4.388 4.388
observation->EnspOLF-1994                  6.645 6.647 6.678 6.678 6.678
observation->EnspOLF-1996                  7.899 7.897 7.906 7.906 7.906
observation->EnspOLF-1998                  5.301 5.301 5.303 5.303 5.303
observation->Enspage-1999                  14.331 14.332 14.287 14.287 14.287
observation->Enspage-2000                  13.973 13.975 13.960 13.960 13.960
observation->Enspage-2001                  13.959 13.979 13.929 13.929 13.929
observation->Enspage-2002                  12.087 12.081 12.078 12.078 12.078
observation->Enspage-2003                  9.702 9.700 9.706 9.706 9.706
observation->Enspage-2004                  12.972 12.966 12.986 12.986 12.986
observation->Enspage-2005                  11.009 11.012 11.009 11.009 11.009
observation->Enspage-2006                  12.744 12.741 12.712 12.712 12.712
observation->Enspage-2007                  11.060 11.061 11.052 11.052 11.052
observation->Enspage-2008                  12.100 12.104 12.101 12.101 12.101
observation->Enspage-2009                  10.055 10.056 10.055 10.055 10.055
observation->Enspage-2010                  17.576 17.581 17.627 17.627 17.627
observation->Enspage-2011                  13.402 13.399 13.412 13.412 13.412
observation->Enspage-2012                  13.148 13.149 13.134 13.134 13.134
observation->Enspage-2013                  16.691 16.694 16.659 16.659 16.659
observation->Enspage-2014                  13.868 13.869 13.870 13.870 13.870
observation->Enspage-2015                  12.764 12.760 12.776 12.776 12.776
observation->Enspage-2016                  19.222 19.224 19.203 19.203 19.203
observation->Enspage-2017                  15.552 15.562 15.511 15.511 15.511
observation->WnspOLF-1992                  10.124 10.122 10.129 10.129 10.129
observation->WnspOLF-1993                  9.047 9.045 8.991 8.991 8.991
observation->WnspOLF-1994                  15.931 15.930 15.893 15.893 15.893
observation->WnspOLF-1996                  9.784 9.788 9.752 9.752 9.752
observation->WnspOLF-1999                  11.553 11.553 11.593 11.593 11.593
observation->WnspOLF-2000                  13.817 13.812 13.892 13.892 13.892
observation->Wnspage-2001                  8.036 8.032 7.994 7.994 7.994
observation->Wnspage-2002                  8.714 8.715 8.701 8.701 8.701
observation->Wnspage-2003                  10.473 10.473 10.507 10.507 10.507
observation->Wnspage-2004                  7.154 7.155 7.140 7.140 7.140
observation->Wnspage-2006                  9.590 9.589 9.551 9.551 9.551
observation->Wnspage-2007                  10.860 10.861 10.822 10.822 10.822
observation->Wnspage-2008                  12.422 12.425 12.435 12.435 12.435
observation->Wnspage-2009                  14.719 14.716 14.790 14.790 14.790
observation->Wnspage-2010                  11.405 11.405 11.370 11.370 11.370
observation->Wnspage-2011                  11.265 11.267 11.275 11.275 11.275
observation->Wnspage-2012                  10.453 10.454 10.465 10.465 10.465
observation->Wnspage-2013                  13.634 13.632 13.634 13.634 13.634
observation->Wnspage-2014                  12.609 12.608 12.597 12.597 12.597
observation->Wnspage-2016                  15.683 15.687 15.722 15.722 15.722
observation->CRsumage-1992                  15.570 15.584 15.566 15.566 15.566
observation->CRsumage-1993                  10.613 10.616 10.624 10.624 10.624
observation->CRsumage-1994                  18.174 18.170 18.130 18.130 18.130
observation->CRsumage-1995                  14.287 14.285 14.242 14.242 14.242
observation->CRsumage-1996                  11.566 11.567 11.607 11.607 11.607
observation->CRsumage-1997                  15.750 15.754 15.688 15.688 15.688
observation->CRsumage-1998                  14.686 14.690 14.690 14.690 14.690
observation->CRsumage-1999                  12.061 12.061 12.084 12.084 12.084
observation->CRsumage-2000                  14.326 14.323 14.283 14.283 14.283
observation->CRsumage-2001                  19.651 19.659 19.520 19.520 19.520
observation->CRsumage-2002                  10.971 10.969 10.960 10.960 10.960
observation->CRsumage-2003                  11.254 11.251 11.210 11.210 11.210
observation->CRsumage-2004                  9.556 9.553 9.561 9.561 9.561
observation->CRsumage-2005                  12.749 12.753 12.654 12.654 12.654
observation->CRsumage-2006                  10.317 10.316 10.317 10.317 10.317
observation->CRsumage-2007                  13.561 13.562 13.504 13.504 13.504
observation->CRsumage-2008                  11.665 11.661 11.639 11.639 11.639
observation->CRsumage-2009                  12.002 12.006 11.965 11.965 11.965
observation->CRsumage-2010                  11.052 11.048 11.055 11.055 11.055
observation->CRsumage-2011                  12.577 12.579 12.571 12.571 12.571
observation->CRsumage-2012                  17.230 17.228 17.141 17.141 17.141
observation->CRsumage-2013                  11.845 11.843 11.822 11.822 11.822
observation->CRsumage-2014                  11.288 11.286 11.271 11.271 11.271
observation->CRsumage-2016                  9.867 9.866 9.847 9.847 9.847
observation->CRsumage-2018                  11.935 11.931 11.895 11.895 11.895
observation->SAsumage-1992                  12.241 12.246 12.269 12.269 12.269
observation->SAsumage-1993                  12.483 12.470 12.502 12.502 12.502
observation->SAsumage-1994                  10.808 10.811 10.818 10.818 10.818
observation->SAsumage-2001                  12.762 12.749 12.839 12.839 12.839
observation->SAsumage-2002                  12.216 12.209 12.254 12.254 12.254
observation->SAsumage-2003                  13.388 13.386 13.386 13.386 13.386
observation->SAsumage-2004                  11.064 11.067 11.005 11.005 11.005
observation->SAsumage-2005                  15.280 15.278 15.386 15.386 15.386
observation->SAsumage-2006                  14.110 14.112 14.078 14.078 14.078
observation->SAsumage-2007                  10.649 10.649 10.639 10.639 10.639
observation->SAsumage-2008                  13.766 13.771 13.836 13.836 13.836
observation->SAsumage-2009                  15.520 15.513 15.530 15.530 15.530
observation->SAsumage-2010                  15.082 15.081 15.130 15.130 15.130
observation->SAsumage-2012                  13.439 13.437 13.479 13.479 13.479
observation->SAsumage-2013                  11.956 11.954 11.982 11.982 11.982
observation->SAsumage-2015                  13.524 13.522 13.540 13.540 13.540
observation->SAsumage-2017                  12.487 12.492 12.478 12.478 12.478
observation->SAautage-1992                  11.819 11.809 11.857 11.857 11.857
observation->SAautage-1996                  10.450 10.449 10.426 10.426 10.426
observation->SAautage-1998                  12.810 12.809 12.835 12.835 12.835
observation->pspawn_1993                  -6.016 -6.017 -6.015 -6.015 -6.015
observation->pspawn-1992                  -4.852 -4.848 -4.829 -4.829 -4.829
observation->pspawn-1998                  -8.291 -8.284 -8.300 -8.300 -8.300
prior->CSacousq->catchability[CSacous].q                -0.847 -0.849 -0.839 -0.839 -0.839
prior->WCacousq->catchability[WCacous].q                -1.163 -1.162 -1.202 -1.202 -1.202
prior->CRsumq->catchability[CRsum].q                -2.153 -2.155 -2.207 -2.207 -2.207
prior->SAsumq->catchability[SAsum].q                -1.760 -1.762 -1.557 -1.557 -1.557
prior->SAautq->catchability[SAaut].q                -1.911 -1.913 -1.762 -1.762 -1.762
prior->CR_process_error->observation[CRsumbio].process_error    0.000 0.000 0.000 0.000 0.000
prior->SA_process_error->observation[SAsumbio].process_error    0.000 0.000 0.000 0.000 0.000
prior->B0_E_with_total_log_b0_prior->process[recruit_E].b0               0.000 0.000 0.000 0.000 0.000
prior->B0_W_with_proportion_prior->process[recruit_W].b0               12.545 12.544 12.555 12.555 12.555
prior->M_all->process[Instant_mortality].m{west.sa}              -1.083 -1.085 -1.217 -1.217 -1.217
prior->Wspsl_shift_param->time_varying[shifted_mu].a                0.422 0.420 0.378 0.378 0.378
prior->YCS_E->process[recruit_E].ycs_values{1975} -0.624 -0.624 -0.623 -0.623 -0.623
prior->YCS_E->process[recruit_E].ycs_values{1976} -0.314 -0.316 -0.319 -0.319 -0.319
prior->YCS_E->process[recruit_E].ycs_values{1977} -0.225 -0.227 -0.238 -0.238 -0.238
prior->YCS_E->process[recruit_E].ycs_values{1978} 0.656 0.654 0.615 0.615 0.615
prior->YCS_E->process[recruit_E].ycs_values{1979} 0.086 0.084 0.062 0.062 0.062
prior->YCS_E->process[recruit_E].ycs_values{1980} -0.630 -0.630 -0.627 -0.627 -0.627
prior->YCS_E->process[recruit_E].ycs_values{1981} -0.619 -0.619 -0.623 -0.623 -0.623
prior->YCS_E->process[recruit_E].ycs_values{1982} -0.551 -0.550 -0.542 -0.542 -0.542
prior->YCS_E->process[recruit_E].ycs_values{1983} -0.627 -0.627 -0.630 -0.630 -0.630
prior->YCS_E->process[recruit_E].ycs_values{1984} -0.591 -0.591 -0.596 -0.596 -0.596
prior->YCS_E->process[recruit_E].ycs_values{1985} -0.641 -0.641 -0.642 -0.642 -0.642
prior->YCS_E->process[recruit_E].ycs_values{1986} -0.408 -0.408 -0.400 -0.400 -0.400
prior->YCS_E->process[recruit_E].ycs_values{1987} 0.142 0.142 0.135 0.135 0.135
prior->YCS_E->process[recruit_E].ycs_values{1988} 0.025 0.026 0.019 0.019 0.019
prior->YCS_E->process[recruit_E].ycs_values{1989} -0.631 -0.631 -0.632 -0.632 -0.632
prior->YCS_E->process[recruit_E].ycs_values{1990} 0.159 0.159 0.158 0.158 0.158
prior->YCS_E->process[recruit_E].ycs_values{1991} -0.643 -0.643 -0.643 -0.643 -0.643
prior->YCS_E->process[recruit_E].ycs_values{1992} -0.643 -0.643 -0.643 -0.643 -0.643
prior->YCS_E->process[recruit_E].ycs_values{1993} -0.643 -0.643 -0.643 -0.643 -0.643
prior->YCS_E->process[recruit_E].ycs_values{1994} -0.531 -0.530 -0.525 -0.525 -0.525
prior->YCS_E->process[recruit_E].ycs_values{1995} -0.534 -0.534 -0.538 -0.538 -0.538
prior->YCS_E->process[recruit_E].ycs_values{1996} -0.573 -0.574 -0.577 -0.577 -0.577
prior->YCS_E->process[recruit_E].ycs_values{1997} -0.623 -0.623 -0.622 -0.622 -0.622
prior->YCS_E->process[recruit_E].ycs_values{1998} -0.511 -0.512 -0.513 -0.513 -0.513
prior->YCS_E->process[recruit_E].ycs_values{1999} 2.014 2.014 2.014 2.014 2.014
prior->YCS_E->process[recruit_E].ycs_values{2000} -0.616 -0.616 -0.614 -0.614 -0.614
prior->YCS_E->process[recruit_E].ycs_values{2001} -0.065 -0.065 -0.071 -0.071 -0.071
prior->YCS_E->process[recruit_E].ycs_values{2002} -0.542 -0.543 -0.545 -0.545 -0.545
prior->YCS_E->process[recruit_E].ycs_values{2003} -0.586 -0.586 -0.585 -0.585 -0.585
prior->YCS_E->process[recruit_E].ycs_values{2004} -0.568 -0.568 -0.568 -0.568 -0.568
prior->YCS_E->process[recruit_E].ycs_values{2005} -0.579 -0.579 -0.578 -0.578 -0.578
prior->YCS_E->process[recruit_E].ycs_values{2006} -0.483 -0.483 -0.482 -0.482 -0.482
prior->YCS_E->process[recruit_E].ycs_values{2007} -0.567 -0.567 -0.566 -0.566 -0.566
prior->YCS_E->process[recruit_E].ycs_values{2008} -0.643 -0.643 -0.643 -0.643 -0.643
prior->YCS_E->process[recruit_E].ycs_values{2009} -0.561 -0.561 -0.562 -0.562 -0.562
prior->YCS_E->process[recruit_E].ycs_values{2010} -0.031 -0.031 -0.046 -0.046 -0.046
prior->YCS_E->process[recruit_E].ycs_values{2011} -0.408 -0.408 -0.408 -0.408 -0.408
prior->YCS_E->process[recruit_E].ycs_values{2012} -0.019 -0.019 -0.023 -0.023 -0.023
prior->YCS_E->process[recruit_E].ycs_values{2013} -0.238 -0.238 -0.242 -0.242 -0.242
prior->YCS_E->process[recruit_E].ycs_values{2014} -0.317 -0.317 -0.314 -0.314 -0.314
prior->YCS_E->process[recruit_E].ycs_values{2015} -0.284 -0.279 -0.302 -0.302 -0.302
prior->YCS_E->process[recruit_E].ycs_values{2016} 0.364 0.363 0.373 0.373 0.373
prior->YCS_W->process[recruit_W].ycs_values{1975} -0.616 -0.616 -0.617 -0.617 -0.617
prior->YCS_W->process[recruit_W].ycs_values{1976} -0.195 -0.195 -0.184 -0.184 -0.184
prior->YCS_W->process[recruit_W].ycs_values{1977} -0.376 -0.376 -0.396 -0.396 -0.396
prior->YCS_W->process[recruit_W].ycs_values{1978} -0.397 -0.396 -0.417 -0.417 -0.417
prior->YCS_W->process[recruit_W].ycs_values{1979} -0.607 -0.607 -0.612 -0.612 -0.612
prior->YCS_W->process[recruit_W].ycs_values{1980} -0.568 -0.567 -0.573 -0.573 -0.573
prior->YCS_W->process[recruit_W].ycs_values{1981} -0.502 -0.502 -0.508 -0.508 -0.508
prior->YCS_W->process[recruit_W].ycs_values{1982} -0.606 -0.606 -0.607 -0.607 -0.607
prior->YCS_W->process[recruit_W].ycs_values{1983} -0.108 -0.109 -0.111 -0.111 -0.111
prior->YCS_W->process[recruit_W].ycs_values{1984} -0.161 -0.163 -0.163 -0.163 -0.163
prior->YCS_W->process[recruit_W].ycs_values{1985} -0.619 -0.619 -0.620 -0.620 -0.620
prior->YCS_W->process[recruit_W].ycs_values{1986} -0.097 -0.098 -0.098 -0.098 -0.098
prior->YCS_W->process[recruit_W].ycs_values{1987} 0.789 0.788 0.782 0.782 0.782
prior->YCS_W->process[recruit_W].ycs_values{1988} -0.251 -0.252 -0.255 -0.255 -0.255
prior->YCS_W->process[recruit_W].ycs_values{1989} 1.207 1.207 1.211 1.211 1.211
prior->YCS_W->process[recruit_W].ycs_values{1990} 1.198 1.197 1.204 1.204 1.204
prior->YCS_W->process[recruit_W].ycs_values{1991} -0.081 -0.081 -0.080 -0.080 -0.080
prior->YCS_W->process[recruit_W].ycs_values{1992} 0.480 0.480 0.480 0.480 0.480
prior->YCS_W->process[recruit_W].ycs_values{1993} -0.643 -0.643 -0.643 -0.643 -0.643
prior->YCS_W->process[recruit_W].ycs_values{1994} 0.014 0.014 0.020 0.020 0.020
prior->YCS_W->process[recruit_W].ycs_values{1995} -0.409 -0.410 -0.413 -0.413 -0.413
prior->YCS_W->process[recruit_W].ycs_values{1996} -0.301 -0.301 -0.306 -0.306 -0.306
prior->YCS_W->process[recruit_W].ycs_values{1997} -0.467 -0.466 -0.471 -0.471 -0.471
prior->YCS_W->process[recruit_W].ycs_values{1998} -0.323 -0.322 -0.334 -0.334 -0.334
prior->YCS_W->process[recruit_W].ycs_values{1999} 2.014 2.014 2.014 2.014 2.014
prior->YCS_W->process[recruit_W].ycs_values{2000} -0.643 -0.643 -0.643 -0.643 -0.643
prior->YCS_W->process[recruit_W].ycs_values{2001} 0.855 0.856 0.846 0.846 0.846
prior->YCS_W->process[recruit_W].ycs_values{2002} -0.222 -0.222 -0.222 -0.222 -0.222
prior->YCS_W->process[recruit_W].ycs_values{2003} -0.612 -0.612 -0.612 -0.612 -0.612
prior->YCS_W->process[recruit_W].ycs_values{2004} -0.570 -0.570 -0.571 -0.571 -0.571
prior->YCS_W->process[recruit_W].ycs_values{2005} -0.574 -0.574 -0.576 -0.576 -0.576
prior->YCS_W->process[recruit_W].ycs_values{2006} -0.280 -0.280 -0.279 -0.279 -0.279
prior->YCS_W->process[recruit_W].ycs_values{2007} -0.438 -0.438 -0.436 -0.436 -0.436
prior->YCS_W->process[recruit_W].ycs_values{2008} -0.464 -0.464 -0.463 -0.463 -0.463
prior->YCS_W->process[recruit_W].ycs_values{2009} -0.514 -0.514 -0.508 -0.508 -0.508
prior->YCS_W->process[recruit_W].ycs_values{2010} 0.070 0.071 0.078 0.078 0.078
prior->YCS_W->process[recruit_W].ycs_values{2011} 0.030 0.031 0.043 0.043 0.043
prior->YCS_W->process[recruit_W].ycs_values{2012} -0.049 -0.050 -0.059 -0.059 -0.059
prior->YCS_W->process[recruit_W].ycs_values{2013} -0.488 -0.489 -0.490 -0.490 -0.490
prior->YCS_W->process[recruit_W].ycs_values{2014} 1.727 1.730 1.742 1.742 1.742
prior->YCS_W->process[recruit_W].ycs_values{2015} 1.684 1.681 1.687 1.687 1.687
prior->YCS_W->process[recruit_W].ycs_values{2016} 0.364 0.363 0.373 0.373 0.373
prior->sel_Whome->selectivity[sel_Whome].v{1}             0.000 0.000 0.000 0.000 0.000
prior->sel_Whome->selectivity[sel_Whome].v{2}             0.000 0.000 0.000 0.000 0.000
prior->sel_Whome->selectivity[sel_Whome].v{3}             0.000 0.000 0.000 0.000 0.000
prior->sel_Whome->selectivity[sel_Whome].v{4}             0.000 0.000 0.000 0.000 0.000
prior->sel_Whome->selectivity[sel_Whome].v{5}             0.000 0.000 0.000 0.000 0.000
prior->sel_Whome->selectivity[sel_Whome].v{6}             0.000 0.000 0.000 0.000 0.000
prior->sel_Whome->selectivity[sel_Whome].v{7}             0.000 0.000 0.000 0.000 0.000
prior->sel_Espmg_male->selectivity[sel_Espmg].v{1}             0.000 0.000 0.000 0.000 0.000
prior->sel_Espmg_male->selectivity[sel_Espmg].v{2}             0.000 0.000 0.000 0.000 0.000
prior->sel_Espmg_male->selectivity[sel_Espmg].v{3}             0.000 0.000 0.000 0.000 0.000
prior->sel_Espmg_male->selectivity[sel_Espmg].v{4}             0.000 0.000 0.000 0.000 0.000
prior->sel_Espmg_male->selectivity[sel_Espmg].v{5}             0.000 0.000 0.000 0.000 0.000
prior->sel_Espmg_male->selectivity[sel_Espmg].v{6}             0.000 0.000 0.000 0.000 0.000
prior->sel_Espmg_male->selectivity[sel_Espmg].v{7}             0.000 0.000 0.000 0.000 0.000
prior->sel_Espmg_male->selectivity[sel_Espmg].v{8}             0.000 0.000 0.000 0.000 0.000
prior->sel_Wspmg_male->selectivity[sel_Wspmg].v{1}             0.000 0.000 0.000 0.000 0.000
prior->sel_Wspmg_male->selectivity[sel_Wspmg].v{2}             0.000 0.000 0.000 0.000 0.000
prior->sel_Wspmg_male->selectivity[sel_Wspmg].v{3}             0.000 0.000 0.000 0.000 0.000
prior->sel_Wspmg_male->selectivity[sel_Wspmg].v{4}             0.000 0.000 0.000 0.000 0.000
prior->sel_Wspmg_male->selectivity[sel_Wspmg].v{5}             0.000 0.000 0.000 0.000 0.000
prior->sel_Wspmg_male->selectivity[sel_Wspmg].v{6}             0.000 0.000 0.000 0.000 0.000
prior->sel_Wspmg_male->selectivity[sel_Wspmg].v{7}             0.000 0.000 0.000 0.000 0.000
prior->sel_Wspmg_male->selectivity[sel_Wspmg].v{8}             0.000 0.000 0.000 0.000 0.000
prior->Enspsl_mu->selectivity[Enspsl].mu               0.000 0.000 0.000 0.000 0.000
prior->Enspsl_s_l->selectivity[Enspsl].sigma_l          0.000 0.000 0.000 0.000 0.000
prior->Enspsl_s_r->selectivity[Enspsl].sigma_r          0.000 0.000 0.000 0.000 0.000
prior->Wnspsl_mu->selectivity[Wnspsl].mu               0.000 0.000 0.000 0.000 0.000
prior->Wnspsl_s_l->selectivity[Wnspsl].sigma_l          0.000 0.000 0.000 0.000 0.000
prior->Wnspsl_s_r->selectivity[Wnspsl].sigma_r          0.000 0.000 0.000 0.000 0.000
prior->Espsl_mu->selectivity[Espsl].mu               0.000 0.000 0.000 0.000 0.000
prior->Espsl_s_l->selectivity[Espsl].sigma_l          0.000 0.000 0.000 0.000 0.000
prior->Espsl_s_r->selectivity[Espsl].sigma_r          0.000 0.000 0.000 0.000 0.000
prior->CRsl_mu->selectivity[CRsl].mu               0.000 0.000 0.000 0.000 0.000
prior->CRsl_s_l->selectivity[CRsl].sigma_l          0.000 0.000 0.000 0.000 0.000
prior->CRsl_s_r->selectivity[CRsl].sigma_r          0.000 0.000 0.000 0.000 0.000
prior->SAsl_mu->selectivity[SAsl].mu               0.000 0.000 0.000 0.000 0.000
prior->SAsl_s_l->selectivity[SAsl].sigma_l          0.000 0.000 0.000 0.000 0.000
prior->SAsl_s_r->selectivity[SAsl].sigma_r          0.000 0.000 0.000 0.000 0.000
additional_prior->YCS.eq.16               0.000 0.000 0.000 0.000 0.000
additional_prior->sp.migr             0.009 0.009 0.013 0.013 0.013
jacobian->log_sum_b0                  0.000 0.000 0.000 0.000 0.000
total_score                  1946.800 1946.810 1946.540 1946.540 1946.540

## [1] "CASAL base model convergence information"
## [1] "fmm: have converged: t = 8.24015e-007 f = 1946.54"                                                        
## [2] "Successful convergence in optimise"                                                                       
## [3] "Minimiser achieved convergence after 599 quasi-Newton iterations using 669 objective function evaluations"
## [1] ""
## [1] "CASAL sensitivity 1 model convergence information"
## [1] "fmm: step size too small. Indicates successful convergence (though this is not the textbook ideal convergence situation)"
## [2] "Successful convergence in optimise"                                                                                      
## [3] "Minimiser achieved convergence after 665 quasi-Newton iterations using 765 objective function evaluations"
## [1] ""
## [1] "Casal2 betadiff_casal_flags_on model convergence information"
## [1] "fmm: have converged: t = 9.38867e-07 f = 1946.8"
## [2] "Successful convergence in optimise"             
## [1] ""
## [1] "Casal2 betadiff_casal_flags_off model convergence information"
## [1] "fmm: have converged: t = 8.80274e-07 f = 1946.81"
## [2] "Successful convergence in optimise"              
## [1] ""
## [1] "Casal2 betadiff_casal_flags_on_low_tol model convergence information"
## [1] "fmm: step size too small. Indicates successful convergence (though this is not the textbook ideal convergence situation)"
## [2] "Successful convergence in optimise"                                                                                      
## [1] ""
## [1] "Casal2 cppad_casal_flags_on model convergence information"
##  [1] "Number of nonzeros in equality constraint Jacobian...:        0"
##  [2] "Number of nonzeros in inequality constraint Jacobian.:      133"
##  [3] "Number of nonzeros in Lagrangian Hessian.............:     8911"
##  [4] "Number of Iterations....: 56"                                   
##  [5] "Number of objective function evaluations             = 92"      
##  [6] "Number of objective gradient evaluations             = 57"      
##  [7] "Number of equality constraint evaluations            = 0"       
##  [8] "Number of inequality constraint evaluations          = 92"      
##  [9] "Number of equality constraint Jacobian evaluations   = 0"       
## [10] "Number of inequality constraint Jacobian evaluations = 57"      
## [11] "Number of Lagrangian Hessian evaluations             = 56"      
## [12] "EXIT: Optimal Solution Found."                                  
## [1] ""
## [1] "Casal2 cppad_casal_flags_off model convergence information"
##  [1] "Number of nonzeros in equality constraint Jacobian...:        0"
##  [2] "Number of nonzeros in inequality constraint Jacobian.:      133"
##  [3] "Number of nonzeros in Lagrangian Hessian.............:     8911"
##  [4] "Number of Iterations....: 56"                                   
##  [5] "Number of objective function evaluations             = 95"      
##  [6] "Number of objective gradient evaluations             = 57"      
##  [7] "Number of equality constraint evaluations            = 0"       
##  [8] "Number of inequality constraint evaluations          = 95"      
##  [9] "Number of equality constraint Jacobian evaluations   = 0"       
## [10] "Number of inequality constraint Jacobian evaluations = 57"      
## [11] "Number of Lagrangian Hessian evaluations             = 56"      
## [12] "EXIT: Optimal Solution Found."
## [1] "CASAL model warnings"
##                  parameter estimate lower.bound upper.bound
## 1   recruitment[E].YCS[25]     0.06        0.06         8.6
## 2   recruitment[W].YCS[25]     0.06        0.06         8.6
## 3 selectivity[CRsl].all[1]        1           1          17
## 4 selectivity[CRsl].all[2]        1           1          50
## [1] ""
## [1] "CASAL sensitivity 1 model warnings"
##                  parameter estimate lower.bound upper.bound
## 1   recruitment[E].YCS[25]     0.06        0.06         8.6
## 2   recruitment[W].YCS[25]     0.06        0.06         8.6
## 3 selectivity[CRsl].all[1]        1           1          17
## 4 selectivity[CRsl].all[2]        1           1          50
## [1] ""
## [1] "Casal2 betadiff_casal_flags_on model warnings"
## $warnings_found
## [1] 8
## 
## $warning_0
## [1] "Estimates were removed because of matching lower and upper bounds. Originally had 134 estimates, now have 133"
## 
## $warning_1
## [1] "estimated parameter 'process[recruit_E].b0' was within 0.001 of upper bound 16.2"
## 
## $warning_2
## [1] "estimated parameter 'process[recruit_W].b0' was within 0.001 of upper bound 0.59"
## 
## $warning_3
## [1] "estimated parameter 'process[recruit_E].ycs_values{1999}' was within 0.001 of lower bound 0.06"
## 
## $warning_4
## [1] "estimated parameter 'process[recruit_W].ycs_values{1999}' was within 0.001 of lower bound 0.06"
## 
## $warning_5
## [1] "estimated parameter 'selectivity[Enspsl].sigma_r' was within 0.001 of upper bound 50"
## 
## $warning_6
## [1] "estimated parameter 'selectivity[CRsl].mu' was within 0.001 of lower bound 1"
## 
## $warning_7
## [1] "estimated parameter 'selectivity[SAsl].sigma_r' was within 0.001 of lower bound 1"
## 
## $type
## [1] "warnings"
## 
## [1] ""
## [1] "Casal2 betadiff_casal_flags_off model warnings"
## $warnings_found
## [1] 10
## 
## $warning_0
## [1] "Estimates were removed because of matching lower and upper bounds. Originally had 134 estimates, now have 133"
## 
## $warning_1
## [1] "estimated parameter 'process[recruit_E].b0' was within 0.001 of upper bound 16.2"
## 
## $warning_2
## [1] "estimated parameter 'process[recruit_W].b0' was within 0.001 of upper bound 0.59"
## 
## $warning_3
## [1] "estimated parameter 'process[recruit_E].ycs_values{1999}' was within 0.001 of lower bound 0.06"
## 
## $warning_4
## [1] "estimated parameter 'process[recruit_W].ycs_values{1999}' was within 0.001 of lower bound 0.06"
## 
## $warning_5
## [1] "estimated parameter 'selectivity[sel_Whome].v{7}' was within 0.001 of upper bound 1"
## 
## $warning_6
## [1] "estimated parameter 'selectivity[Enspsl].sigma_r' was within 0.001 of upper bound 50"
## 
## $warning_7
## [1] "estimated parameter 'selectivity[Wnspsl].sigma_r' was within 0.001 of upper bound 50"
## 
## $warning_8
## [1] "estimated parameter 'selectivity[CRsl].mu' was within 0.001 of lower bound 1"
## 
## $warning_9
## [1] "estimated parameter 'selectivity[SAsl].sigma_r' was within 0.001 of lower bound 1"
## 
## $type
## [1] "warnings"
## 
## [1] ""
## [1] "Casal2 betadiff_casal_flags_on_low_tol model warnings"
## $warnings_found
## [1] 6
## 
## $warning_0
## [1] "Estimates were removed because of matching lower and upper bounds. Originally had 134 estimates, now have 133"
## 
## $warning_1
## [1] "estimated parameter 'process[recruit_E].b0' was within 0.001 of upper bound 16.2"
## 
## $warning_2
## [1] "estimated parameter 'process[recruit_W].b0' was within 0.001 of upper bound 0.59"
## 
## $warning_3
## [1] "estimated parameter 'process[recruit_E].ycs_values{1999}' was within 0.001 of lower bound 0.06"
## 
## $warning_4
## [1] "estimated parameter 'process[recruit_W].ycs_values{1999}' was within 0.001 of lower bound 0.06"
## 
## $warning_5
## [1] "estimated parameter 'selectivity[CRsl].mu' was within 0.001 of lower bound 1"
## 
## $type
## [1] "warnings"
## 
## [1] ""
## [1] "Casal2 cppad_casal_flags_on model warnings"
## $warnings_found
## [1] 6
## 
## $warning_0
## [1] "Estimates were removed because of matching lower and upper bounds. Originally had 134 estimates, now have 133"
## 
## $warning_1
## [1] "estimated parameter 'process[recruit_E].b0' was within 0.001 of upper bound 16.2"
## 
## $warning_2
## [1] "estimated parameter 'process[recruit_W].b0' was within 0.001 of upper bound 0.59"
## 
## $warning_3
## [1] "estimated parameter 'process[recruit_E].ycs_values{1999}' was within 0.001 of lower bound 0.06"
## 
## $warning_4
## [1] "estimated parameter 'process[recruit_W].ycs_values{1999}' was within 0.001 of lower bound 0.06"
## 
## $warning_5
## [1] "estimated parameter 'selectivity[CRsl].mu' was within 0.001 of lower bound 1"
## 
## $type
## [1] "warnings"
## 
## [1] ""
## [1] "Casal2 cppad_casal_flags_off model warnings"
## $warnings_found
## [1] 6
## 
## $warning_0
## [1] "Estimates were removed because of matching lower and upper bounds. Originally had 134 estimates, now have 133"
## 
## $warning_1
## [1] "estimated parameter 'process[recruit_E].b0' was within 0.001 of upper bound 16.2"
## 
## $warning_2
## [1] "estimated parameter 'process[recruit_W].b0' was within 0.001 of upper bound 0.59"
## 
## $warning_3
## [1] "estimated parameter 'process[recruit_E].ycs_values{1999}' was within 0.001 of lower bound 0.06"
## 
## $warning_4
## [1] "estimated parameter 'process[recruit_W].ycs_values{1999}' was within 0.001 of lower bound 0.06"
## 
## $warning_5
## [1] "estimated parameter 'selectivity[CRsl].mu' was within 0.001 of lower bound 1"
## 
## $type
## [1] "warnings"

1.4 Matching of outputs

Time series comparisons with CASAL base model results

## [1] "Catch time series base model comparison for run betadiff_casal_flags_on"
## [1] "Actual catches for Ensp1 match: yes"
## [1] "Actual catches for Wnsp1 match: yes"
## [1] "Actual catches for Ensp2 match: yes"
## [1] "Actual catches for Wnsp2 match: yes"
## [1] "Actual catches for Esp match: yes"
## [1] "Actual catches for Wsp match: yes"
## [1] ""
## [1] "Catch time series base model comparison for run betadiff_casal_flags_off"
## [1] "Actual catches for Ensp1 match: yes"
## [1] "Actual catches for Wnsp1 match: yes"
## [1] "Actual catches for Ensp2 match: yes"
## [1] "Actual catches for Wnsp2 match: yes"
## [1] "Actual catches for Esp match: yes"
## [1] "Actual catches for Wsp match: yes"
## [1] ""
## [1] "Catch time series base model comparison for run betadiff_casal_flags_on_low_tol"
## [1] "Actual catches for Ensp1 match: yes"
## [1] "Actual catches for Wnsp1 match: yes"
## [1] "Actual catches for Ensp2 match: yes"
## [1] "Actual catches for Wnsp2 match: yes"
## [1] "Actual catches for Esp match: yes"
## [1] "Actual catches for Wsp match: yes"
## [1] ""
## [1] "Catch time series base model comparison for run cppad_casal_flags_on"
## [1] "Actual catches for Ensp1 match: yes"
## [1] "Actual catches for Wnsp1 match: yes"
## [1] "Actual catches for Ensp2 match: yes"
## [1] "Actual catches for Wnsp2 match: yes"
## [1] "Actual catches for Esp match: yes"
## [1] "Actual catches for Wsp match: yes"
## [1] ""
## [1] "Catch time series base model comparison for run cppad_casal_flags_off"
## [1] "Actual catches for Ensp1 match: yes"
## [1] "Actual catches for Wnsp1 match: yes"
## [1] "Actual catches for Ensp2 match: yes"
## [1] "Actual catches for Wnsp2 match: yes"
## [1] "Actual catches for Esp match: yes"
## [1] "Actual catches for Wsp match: yes"
## [1] ""

Derived quantities

SB0, SBcurrent, MSY, F_MSY, others…

1.5 Plots

Comparison plots

## [1] "CASAL base parameter correlation range (excluding 1.0): -0.991373 0.996882"

## [1] "CASAL sensitivity 1 parameter correlation range (excluding 1.0): -0.932422 0.996848"

## [1] "Casal2 betadiff_casal_flags_on parameter correlation range (excluding 1.0): -0.889617003107542 0.996691861929604"

## [1] "Casal2 betadiff_casal_flags_off parameter correlation range (excluding 1.0): -0.889151569029662 0.996800894132685"

## [1] "Casal2 betadiff_casal_flags_on_low_tol parameter correlation range (excluding 1.0): -0.921980407455868 0.996701916513473"

## [1] "Casal2 cppad_casal_flags_on parameter correlation range (excluding 1.0): 0 0"

## [1] "Casal2 cppad_casal_flags_off parameter correlation range (excluding 1.0): 0 0"